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Abstract 



Various light-quark channel current-current correlators are subjected to 
the concept of a non-perturbative component of coarse graining in opera- 
Oh, tor product expansions introduced in a parallel work. This procedure allows 

(— | \ for low-energy structure of the OPE-derived spectral function. With naive 

vacuum saturation for 4 quark operators and using lattice data for the gauge 
invariant scalar quark correlator the results are far off the experimentally mea- 
sured behavior. However, using the correlation length of the gauge invariant 
vector quark correlator, which is about 10 times smaller than the scalar one, 
the qualitative results are rather realistic. Namely, the input of information 
on the mass of the lowest resonance in one channel yields the correspond- 
ing masses within acceptable errors in other channels. Still, the shapes of 
the calculated spectral functions are considerably deformed as compared to 
experiment. This may be a consequence of vacuum saturation and the trunca- 
tion at a mass dimension which is below the critical dimension from which on 
the asymptotic expansion does not approximate anymore. To improve on this 
high-resolution lattice information on gauge invariant n > 2 point correlators 
would be needed. Motivated by the small effective correlation length in the 4 
quark contributions the relevance of the approach for heavy quark physics, in 
particular in the calculation of nonleptonic, inclusive Ar, is discussed. 



1 Introduction 



The basic assumption made in QCD sum rules (QSR) in particular [|J and in appli- 
cations of the operator product expansion (OPE) to hadronic physics in general f|| 
is quark-hadron duality, namely the possibility to express the low-energy dispersive 
part of a correlator of QCD currents in terms of hadronic cross sections. For exter- 
nal momentum q with q 2 = —Q 2 < correlators are believed to be (asymptotically) 
approximated by the OPE which is a power series in Q~ 2 . A first estimate for the 
critical dimension d c ~ 12 from which on the entire coefficients, being products of 
Wilson coefficients and averages over local, gauge invariant operators, are dominated 
by short distance effects was performed in ref. PJ for pure gluodynamics. Thereby, 
the vacuum was described by a dilute instanton gas. Powers with d > d c have no 
potential to reflect the hadronic properties of the respective channel which relate to 
the operator averages via dispersion relations in QSR (global duality) and via the 
optical theorem and analytical continuation in applications such as calculations of 
nonleptonic decay widths of heavy mesons (local duality). 

Apart from some puzzlesQ the tremendous success of the sum rule method seems 
to support the global version of quark-hadron duality. 

The observation that an analytical continuation of the OPE to time-like, external 
momenta does not yield the phenomenological resonance structure in the imaginary 
part has lead to the notion of local duality violation 0, [|, [|. The up-to-date 
claim is that the asymptotic nature of the expansion does not allow for additive, 
exponential-like terms to appear. These terms are to be responsible for a "wiggling" 
of the low-energy part of the spectral function . This, however, seems to contradict 
the QSR philosophy [[[]]: If only the first few terms of the OPE (up to the critical 
dimension d c ) are needed to describe the lowest resonances in a given channel via a 
dispersion relation then operators of dimension d > d c should also be forgotten in 
an OPE-based construction of the spectral function. 

Therefore, we propose an alternative in ||. The observation is that the non- 
perturbative behavior of the QCD vacuum, when probed with low euclidean mo- 
menta, is not only characterized by nonvanishing vacuum expectation values (VEV's) 
of local, gauge invariant operators but also by finite correlation lengths of the corre- 
sponding gauge invariant correlators |7|, ||. As a consequence the fundamental field 
operators must loose their relevance with decreasing resolution (see next section). 
From the knowledge of gauge invariant correlation functions of fundamental opera- 
tors at high resolution Q evolution equations for the VEV's A(Q) of the relevant 
local, effective operators can be derived for Q < Q . For d = 4, where 2-point 
functions are needed, we have A(Q) ~ exp[4/5(AQ) _1 ] with A being the correlation 
length. Compared with the conventional power correction the effect becomes noti- 

1 large channel-to-channel variation of spectral continuum thresholds, see ref. ||, large error in 
vacuum averages of local operators, unclarified status of such approximations as vacuum saturation 
in vacuum averages over 4-quark operators, very large non-perturbative effects in the scalar and 
pseudoscalar channels 
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cable if Q ~ A -1 . For gluon and quark condensates the lattice implies A(Q ) y s at 
Qo ~ 2 GeV which are compatible with the phenomenological values obtained from 
conventional QSR's at /i ~ 1 GeV. Since the inverse correlation lengths are below 
the mass of the p-resonance (assuming naive vacuum saturation at dimension 6) 
light-quark channel QCS sum rules are hardly touched by the exponential decrease 
of the operator VEV's ||. However, calculating the spectral function from the trun- 
cated OPE these large correlation lengths do not reflect the resonance physics in 
the respective channels. A much smaller, effective correlation length is needed at 
dimension 6 to yield more realistic spectra. To investigate the properties of OPE- 
based spectral functions in various light-quark channels is the main purpose of this 
paper. 

The presentation is set up as follows: In the next section we briefly review the 
idea of non-perturbative coarse graining of local operator VEV's as it is developed in 
||. Thereby, the focus is on correlation functions which are represented by so-called 
connected diagrams. Section 3 investigates the vacuum saturation hypothesis for 4- 
quark operators |l| in the light of non-perturbative coarse graining. The low-energy 
parts of the spectral functions in the p, ai, it, and (f> channels are calculated in section 
4 by analytical continuation of the OPE to time-like external momenta. Implications 
for the calculation of the difference of nonleptonic inclusive decay widths in neutral 
5-meson systems are discussed in section 5. The last section summarizes the results. 

2 Euclidean exponentials and Minkowskian oscil- 
lations 

In this section we briefly review the work of || on non-perturbative coarse graining 
for VEV's of local, gauge invariant operators. 

At a large euclidean momentum Q = y/Cp, where we expect the description of 
the dynamics in terms of the continuum action and local operators made of funda- 
mental fields to be sufficiently accurate, we start by expanding the current- current 
correlator into a conventional OPE. The evolution to lower momenta is obtained 
by running the Wilson coefficients perturbatively via the running coupling a s (Q 2 ) 
and the anomalous operator dimensions 0. According to || the non-perturbative 
running of an operator average is governed by the non-perturbative part of the corre- 
sponding gauge invariant correlator in euclidean position space. This correspondence 
can be expressed as 
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In eq. parallel transporters 




(2) 
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are appropriately contained in the non-local expression to define gauge invariant 
correlations. The sum runs over all relevant^], piecewise straight f?| trajectories of 
parallel transport, and the symbol V demands path ordering. With the normal- 
ization factor 1/Af, which depends on n and the numbers of fields transforming 
under the fundamental and adjoint representation, the correlation function reduces 
to the "condensate" in the limit x±, - ■ ■ ,x n — > 0. Making the convention that an 
arrow pointing from Xi towards Xj stands for the parallel transport S(xi,Xj), we 
have a way to depict correlators. Note that points with fields transforming under 
the fundamental (adjoint) representation are connected to one (two) lines of parallel 
transport. There are disconnected and connected diagrams. In this section we only 
consider the latter. 

To coarse grain from resolution Q to resolution Q — dQ we average the non- 
perturbative part of the correlator corresponding to a connected diagrams over a 
(euclidean) ball of radius cIRq with 



, p 1 1 1 ( 1 A dQ 



Using the short-hand notation of eq. ([l]), this is written as 

< F <°»«-« = (W)'-» L- m ■ ■ ■ A - {Fii0) ■ ■ ■ F " w> « • (4) 



where 




V(dR Q ) = ^(dR Q y = -^ j^i) . ( r,) 



In Fig. 1 there are 3 examples for connected diagrams. 

Let us now focus on 2-point functions as they are relevant for dimension 3 and 
4 quark and gluon operators, respectively. Only the gauge invariant bilocal quark 
correlator [§] 

(trq{x)S(x,0)q(0)) (6) 
and the gluonic field strength correlator 0, || 

trF^(x)S(x,0)F KX S\x,0)) (7) 

have been measured on the lattice. The results imply that there exists an additive 
decomposition into a perturbative, power-like in \x\, and a non-perturbative, expo- 
nential in \x\ piece 0. We are interested in the non-perturbative part, which in 
short-hand reads 

(F 1 (0)F 2 (x))^ = A(Q)exp(-\x\/\). (8) 

As explained in A is expected to have direct phenomological meaning, and so 
it can not depend on the resolution Q. Hence, coarse graining may only affect the 

2 This will be specified later. 
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x 3 x 2 



(Hi) 

Figure 1: Diagrammatic representation of the gauge invariant correlators corre- 
sponding to the local operators (i) q(0)q(0), (ii) trF fJiU (0)F UK (0)F tt ^(0), and (in) 
tiF^(0)F UK (0)F KX (0)F x ^(0). 
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pre-exponential factor A(Q). An evolution equation for A(Q) can be derived if we 
combine eqs. (f|) and (B): 



-i 



(v(dQ/Q 2 )) " I dQ d*xA(Q)e 



\x\/X 



A(Q - dQ)e 



-0/A 



A(Q - dQ) . (9) 



Note that for A — > oo A is invariant under coarse graining. Therefore, the con- 
ventional treatment of non-perturbative corrections in the framework of the OPE 
ij is recovered, and condensates do not depend on the resolution in this limit. In 
particular, they could be determined at a large resolution, where fundamental fields 
make sense. In the real world, however, correlation lengths are finite |J. We will 
see later how this is reflected in the existence of hadronic resonances. 

Expanding the l.h.s. and r.h.s. of eq. (0) in ^ and comparing coefficients of the 
linear terms, we obtain the following equation: 



^( Q) = -L Q -M(Q), 



with solution 



A(Q) = A(Q )ex V 
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Qo). 



(10) 



11^ 



At dimension 4 there are no anomalous operator dimensions and being content with 
a determination of Wilson coefficients at the lowest possible order in a s the generic 
form of a non-perturbative correction is 



A,(Qo) 
(Q 2 ) 2 



cxp 



1 



5A4 \Q Qq 



(12) 



Eq. ( |TTD implies that compared to conventional dimension 4 power corrections there 
is a noticable suppression if Q is of the order of A -1 or less. 
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Lattice measurements with Np = 4 staggered fermions of mass a ■ m q = 0.01 
and lattice resolution Q = a -1 ~ 2 GeV suggest that scalar fermionic and gluonic 
correlation lengths are A* ~ 3.2 GeV -1 and X g ~ 1.7 GeV -1 , respectively [jSJ. For 
the vector fermionic correlation length A^ ~ 1/10 A* was found ||. So for dimension 
4 and, assuming naive vacuum saturation, also for dimension 6 this does not seem to 
pose a problem for the sum rule analysis of light quark correlators (at Qo = a -1 ~ 2 
GeV A q {Q ) as well as A g (Q ) are compatible with their QSR determined values 
at ii ~ 1 GeV ||). However, as we will see later, naive vacuum saturation gives 
light-quark channel spectra which are completely off the experimentally measured 
behavior. Qualitatively more realistic spectral functions can be obtained using the 
much smaller A^. 

The contribution of dimension 4 to the spectral function p(s) is (up to a nor- 



malization!) obtained by analytically continuing to Q 



[s + ie) or Q = —iy/s, 



(s > 0), and taking the imaginary part. A term like the one in eq. (|T^) corresponds 
to a term 



■ exp 



4 1 
5A4 Qo 



sm 



(13) 



in p(s). The oscillatory behavior manifests itself in a quite different way than it was 
suggested in refs. [§: sin[-y/i 1 ] instead of sin[y / s] or sin[s]. There, oscillations are 
present everywhere though power suppressed at high s. Here, oscillations only start 
if -y/i is larger than the correlation length of the corresponding 2-point function. 



3 Vacuum saturation in the context of gauge in- 
variant correlations 

After the treatment of connected diagrams in the previous section we turn to discon- 
nected diagrams in this section. A sufficient condition for factorized coarse graining 
is the factorization of the non-perturbative part in the corresponding correlation 
function. It is unlikely that such a factorization occurs for connected diagrams. 

Let us focus on the case of 4 quark operators since generalizations to higher 
dimensions are straightforward. Disconnected diagrams are associated with 4 quark 
operators g(0)rg(0)g(0)rg(0) composed of color singlet currents (Fig. 2). Thereby, q 
is a single flavor quark field, and T denotes one of 1, 75, 7 M , • • • or a product of them. 
Diagrams with permutated arguments are irrelevant because of the integrations in 
eq. (f|) and the fact that we may assume translational, P, and T invariance 0. 
In practical applications one encounters color octet currents in 4 quark operators. 
These structures lead to special forms of disconnected diagrams which are due to 
the fact that gauge invariant correlation functions can only be defined if 2 of the 

4 arguments coincide, hence yielding 3 point functions. For example, an operator 
q(0)Tt a q(0)q(0)Tt a q(0) demands non-local contributions of the form 

(t) (q(0)^S(0,x 1 )t a q(x 1 )q(x 2 )S(x 2 ,0)Tt a q(0)) ; 



5 



x 3 x 2 





x l 

Figure 2: Diagrammatic presentation of the gauge invariant correlator corresponding 
to the local operator q(0)Tq(0)q(0)Tq(0). 




*2 *2 



(i) (ii) 

Figure 3: Diagrammatic presentation of the contributions of eq. flUj ) to the gauge 
invariant correlator associated with the local operator g(0)rt a g(0)g(0)rt a g(0). 

(ii) (q(x l )rt a S(x l ,0)q(0)q(x 1 )Ft a S(x 1 ,x 2 )q(x 2 )) . (14) 

Thereby, the color matrices are normalized as tr t a t b = 25 ab . The corresponding 
diagrams are depicted in Fig 3. Lines do not meet at the points (for (i)) and x± 
(for (ii)) because they connect to different fields. 

Since gauge invariant n > 2 point functions have yet not been measured on the 
lattice we do have to think of approximations involving only 2-point functions. The 
(quite intuitive) hypothesis of vacuum saturation can be implemented in two ways: 
1) vacuum saturation on the level of local operators according to the formula of 
ref. [|IJ with separate coarse graining for each factor and 2) vacuum saturation on 
the level of the correlation function, that is, first derealization of the operator and 
then vacuum insertion. There is one obvious reason why the two prescriptions lead to 
different results. In case 1) we essentially square the result for the quark condensate, 
which implies correlations between points of maximal separation dQ/Q 2 , whereas in 
case 2) one factor may contain correlations between points of maximal separation 
2dQ/Q 2 (Fig. 2 or case (ii) in Fig. 3). Another, not-so-obvious reason becomes 
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apparent if we insert the vacuum state into the correlator for singlet currents 

(q(0)TS(0,x 1 )q(x 1 )q(x 2 )TS(x 2 ,x 3 )q(x 3 )) . (15) 

Then information of the 2-point function 

(q{0)rS{0,x)q(x)) (16) 

is needed. From ref.0 we know that on a lattice the correlation length of the 
vector correlator (r = 7^) is about 10 times smaller than the correlation length X s 
of the scalar correlator (r = 1)! This result has been obtained for a quark mass m 
with a-m = 0.01 and a" 1 ~ 2 GeV in the Np = 4 theory with staggered fermions. In 
the octet case vacuum insertion inbetween the currents factorizes the correlator into 
gauge variant 2 point functions. The gauge invariant product of them has not yet 
been measured on the lattice so we can not compare the corresponding correlation 
length with the singlet case. 

Let us keep all these points in mind in the next section, where for simplicity we 
apply prescription 1). With the necessary lattice information available we hope to 
investigate case 2) in a separate publication. 

4 Light-quark channels 

After writing down the coarse grained OPE's in the p, at, tt, and channels, in 
this section we calculate the respective spectral functions by analytical continuation 
from euclidean to time-like momenta. We restrict ourselves to light-quark channels 
because here the two needed correlation functions for corrections up to dimension 6 
(assuming vacuum saturation in the sense of 1)) have been measured on the lattice. 
We investigate with two sets of parameters. Set (A) takes the Np = 4 results of 
refs. H for the lowest quark mass (a ■ m = 0.01) literally, that is 

(A): 

X q = 3.1 GeV" 1 , A^a- 1 = Q ~ 2 GeV) = (0.212 GeV) 3 ; 

X g = 1.7 GeV" 1 , Agia- 1 = Q ~ 2 GeV) = 0.015 (GeV) 4 . (17) 

Due to the uncertainty in the way dimension 6 contributions are treated and mo- 
tivated by the 10 times smaller vector correlation length X v we use X q = X v in 
(B) 

OB): 

X q = 0.3 GeV" 1 , A^ 1 = Q ~ 2 GeV) = (0.212 GeV) 3 ; 

X g = 1.7 GeV" 1 , Agfa- 1 = Q ~ 2 GeV) = 0.015 (GeV) 4 . (18) 
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4.1 /^-correlator 

Here, we investigate the correlator of currents j p = 1/2(m7 m m — d , -f fl d) which at 
Q ~ Qo and up to dimension 6 has the following conventional OPE [l[ 



i f (TtyWjftO)}) = (q»qu-g»uq 2 ) X 



1 + 



a s (Q 2 



57T 



7T 



log Ql + Q' 



- ( vn u uu + m d dd)„ + — ( — F* F a \ 



1 / ' I 



Qo 



\ q=u,d,s I 



After application of the vacuum saturation hypothesis 0, going to the SU(2) chiral 
limit m u = = 0, and implementing the non-perturbative coarse graining of 
operator averages^ the restriction Q ~ Q can be dropped, and the scalar part of 
the correlator (curly brackets on the r.h.s. of eq. (|ll?D) reduces to 




Qo 
(20) 



4.2 ai-correlator 

The vacuum averaged OPE for the correlator of the current j® 1 = 1/2 (^7^75 u — 
o?7/i75 d) at Q ~ Qo and up to dimension 6 reads in the SU(2) chiral limit [|J 

i / rf 4 e^' (T{#On)#(0)}) = - M 2 ) x {- JL (l + ^f^j log|j+ 



O'.s- 



24Q 4 \ vr" ^/ Qo 2Q 



(21) 



3 At dimension 6 there is logarithmic running of the Wilson coefficients due to the nonvanishing 
anomalous dimensions of the corresponding operators. Together with the running coupling a s (Q 2 ) 
these dependences almost cancel [jj], and hence we will omit them throughout what follows. Op- 
erators of dimension 4 are perturbatively invariant under the renormalization group. We take the 
value a s {Ql = (2 GeV) 2 ) - 0.2. This is motivated by ref. |[J where a non-perturbative evolution 
of a s has been obtained in two flavor lattice QCD using the Schrodinger functional scheme. With 
the conversion A^g = 2.382A and taking A^g = 0.4 GeV, one obtains a s (2GeV) 2 ) ~ 0.19. We 
have varied the coupling up to a s = 0.4 but no drastic changes occur. 



8 



After vacuum saturation and implementation of non-perturbative operator running 
the restriction Q ~ Q can be forgotten, and we have 



n ai (Q 2 ) 

na s (Q 2 ) 176 



1 (, , a s (Q 2 ) \i Q 2 , 1 r 



4 , 1 



57T 



exp[ 



n J -Ql 24Q 4 

1 ,_ x 2 



(- - —)}(—F a F a ) 

5A/Q Qo ,l \ir "V 



Q 6 81 r 5A/Q Q ' 
4.3 7r-correlator 

The pion current is given as j* = 1/2 i (1Z75 it — c/75 d), and in the SU(2) chiral limit 
its correlator can be expanded as M 



Qo 
(22) 



i / d^ (T{f{x)f{0)}) = -3Q 5 



1 



16vr 2 



7T i 8 Q§ 48Q 4 \ tt 



Q 6 



\ q=u,d,s I 



After vacuum saturation and with non-perturbative operator running (no restriction 
Q ~ Qo anymore) the piece in curly brackets becomes 



rr(Q 2 ^ 



1 ^ , a s (Q 2 ) \, Q 2 1 r 



16tt 2 V ' vr J "° Qo 48Q 
Q* 81 6XP[ 5A/Q g jJ(gg)Qo ' 



^£ Z7" 17a \ 
5X g y Q Qo'' \ 7T • ' / Qo 



(24) 



4.4 (^-correlator 

The 0-meson current is defined as j^ 1 = — l/3s7^s. If we keep the s-quark mass 
finite, let m u = m d = 0, assume that the s-quark condensate at Q = Q does not 
deviate from the ones of it or d quarks, then at Q ~ Qo the OPE of the ^-correlator 



becomes 1 1 



1 J (T{#(aO#(0)}) = (wh, - g^q 2 ) x 

1 1 /a, 

— (m~nn\ - 4- 

v |_ 

ixa s (Q 2 



2 / W egg 2 

9 I 8vr 2 I 7T 



1 #\ 

los Ql + 



^^0 + 24(7^^ 



Q 6 



(s7 M 7 5 t a s) 2 ) Qo + -(s7 M t a s X! 

\ q=u,d,s I 



(25) 
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m (iiB) (iiiB) (ivB) 

Figure 4: The spectral functions ImU(Q 2 — s — ie), (s > 0) for sets (A) and (B), 
where (i), (ii), (Hi), and (iv) correspond to the p, ai, 7r, and (f> channels, respectively. 



Building in vacuum saturation and non-perturbative running of the operator VEV's, 
the restriction Q ~ Q does not apply anymore, and the scalar piece in curly brackets 
becomes 



exp[— U - — )] (m s qq) Qo + ^exp[- — (- - — )] (^F^ 



1 

™ s (Ql) 112 




4,1 1 „ , , 1 r 4 1 1 v , /a 



5A/Q Qo J yo 24 5A/Q Q \ vr 



exp[ 



g 6 81 5A 



: ^ - ik )] <?9> -) • (26) 



4.5 Spectral functions 

Here, we calculate the spectral functions of the respective channels from the OPE 
by taking the imaginary part of H(Q 2 ) at analytical continued Q 2 = — s — ie or Q = 
\fQ 1 = —iyfs with real s and s > 0. As in the conventional sum rule approach we 
expect to obtain information on the lowest resonances (but now without the apriori 
assumption of quark-hadron duality). Fig. 4 shows the results of the calculations 
for a s (Ql = (2GeV) 2 ) = 0.2 we use m s = 120 MeV. Exact vacuum saturation in 
the sense of 1) is assumed with the commonly introduced "correction" factor k > 1 
set equal to unity. Due to the truncation of the OPE at dimension 6 and other 
errors (large quark mass and Np = 4 in the lattice calculation, vacuum saturation, 
chiral SU(2) limit, and a relatively low "fundamentality scale" Q ) we have to cut 
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off the spectra at some lower bounds. Since spectral functions ought to be positive 
definite natural candidates y/s^ for these cutoffs are the first zeros encountered when 
decreasing yfs. 

Looking at the spectra belonging to set (A) above these cutoffs, there is no 
resemblance of the measured behavior. In contrast, set (B) seems to give a more 
realistic behavior: There is higher spectral strength at lower energy in the p than 
there are in the a\ or <p channels. In the latter we have varied the strange quark 
mass: For m s below 50 MeV the channel behaves p-like, that is, there is a large 
peak at the lower bound of the spectrum. For m s > 50 MeV the 0-spectrum is ax- 
like, that is, the spectral strength is steadily decreasing from its perturbative value 
down to its first zero with decreasing energy. Similarily, the 7r channel exhibits a 
large concentration of spectral strength at low energies. 

Let us process the spectral information contained in Fig. 4 to more quantitative 
statements. One may ask where the "center-of-mass" 



r£(«VS)Im(») 



of a given spectrum within a low-energy domain ^/sq < ^fs < y/s[ is. This is mo- 
tivated by the experimental fact that the lowest resonance dominates the spectrum 
within a large range of energies. Note that using a narrow resonance model, the 
residue of this single resonance cancels in the definition of M (ratio of moments). 
We may now define the universal upper bound y/si of this low-energy domain by 
demanding M p to coincide with the mass m p = 770 MeV of the lowest p-resonance. 

With Jsl = 420 MeV eq. fl27|) then yields ^ = 1.2 GeV. Using this, we calculate 
M for the other channels with always being the first zero of the respective spec- 
tral function. The result can be compared with the measured mass of the lowest 
resonance. We have 

M 855 

= 540 MeV , M ai = 855 MeV , — ^ = — — ~ 0.68 ; 



■ 420 MeV, M w = (>S() MoV . — — — — — : — - (). ( ) I 



1260 

M w _ 680 

l/2[m^ + m^(1300)] ~ 720 



, u = 680 MeV , Ms = 960 MeV , — * = ~ 0.94 . (28) 

ms 1020 

Due to the distinguished role of the 7r-meson as a Goldstone-boson we have taken 
along the next pion resonance (7r(1300)) assuming their residues to be equal. The 
differences between the measured (still using the narrow resonance model) and the 
computed ratios of moments is about 30% for the a% channel and less than 10% 
level for ir and channels. We do emphasize at this point that conventional OPE's 
would have given the same value for M in all channels (apart from small perturbative 
corrections to the Wilson coefficients). 
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Why do the shapes deviate considerably from the experimental ones (as a gen- 
eral feature, they seem to be shifted to lower energies)? Although the use of the 
small correlation length in set (B) was motivated by ambiguities concerning vacuum 
saturation we do believe that the truncation of the OPE at dimension 6 is the major 
source of deviation. It is quite plausible that higher operator dimensions introduce 
shorter and shorter effective correlation lengths and therefore higher mass scales to 
govern the operator VEV's at low resolution ||. Viewed in this context, the use of 
a small correlation length at dimension 6 may be an effective way to simulate higher 
mass dimensions. Apart from this there are, of course, the unresolved problems 
linked to vacuum saturation as such and the way it is being implemented as they 
were discussed in the previous section. 



5 Implications for the calculation of AFb 

Motivated by the occurrence of a large (effective) mass scale in the previous section 
we discuss in this section how the notion of non-perturbative coarse graining of 
local operator averages may influence the calculation of nonleptonic, inclusive width 
differences in neutral 5-meson systems. 

To be specific, we take the example of the B s -B s system as it was treated in 
refs. ||12||. On the level of an effective weak Hamiltonian, which is obtained by in- 



tegrating out the heavy bosons Z°, W of the fundamental theory (also including 
pQCD corrections) and which, omitting Cabibbo suppressed contributions, is of the 
form 



Heff = ^ nVcs (§ Cl ° l + C8 ° S ) ' (29) 



the width difference ATb s between the mass eigenstates 



by virtue of the optical theorem as (relativistic normalization of states) 



Bh/l) can be calculated 



Imi J d 4 xT 7i e f f{x)7i e f f 



B s ) . (30) 



Ar B = — I b s 

s 2M Ba \ 

The first six AB = 1 operators of eq. fl29|) are of the four quark type, for example 
Oi = (b a Cp) V - A (cps a ) V -A, and 8 = g/(8n 2 ) m b b a a^ u {l - lh)F£ v t a aj3 sp. 

To be able to handle eq. (|30"D one expands in position space the T-product into 
an OPE which contains AB = 2 operators, for example 

Q = {baSa)v-A{bpSp) V -A ■ (31) 

In view of the subsequent x integration the OPE only makes sense if the momentum 
q of the external states is in the euclidean (or space-like) domain. So for a meaningful 
calculation formally an analytical continuation of the time-like momentum squared 
q 2 ~ ml of the heavy mesons at rest to q 2 = —Q 2 < must be performed. After 
a determination of the Wilson coefficients in the euclidean the result is analytically 
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continued back to q ~ m . Thereby, one has to deal with two renormahzation scales 



Hi y 2- The former is due to the expansion of the fundamental, electroweak (together 
with QCD corrections) interaction into the effective Hamiltonian, whereas the latter 
comes in via the separate scale dependence of the AB = 2 operators and Wilson 
coefficients. It is claimed in refs. |T2|] that in the AB = 2 OPE the dependence 
on \i\ is almost cancelled on the level of the Wilson coefficients at the same order 



in a s . Since the MS-scheme was used in [T_2| the MS matrix elements had to be 
matched to the matrix elements obtained in lattice renormahzation at the low lattice 
renormahzation point 2~ GeV. Subsequently, they were run up to ~ m&. 

Keeping the scale of the matrix element fixed in a purely perturbative renormal- 
ization group evolution, the transcendental dependences on the external momentum 
scale are powers of logarithms due to their resummation. Violating local duality, 
we have seen in the previous sections that these logarithmic dependences do not 
introduce resonance structure in the spectral functions of light-quark channel vac- 
uum correlators. If the scale of the process is comparable to the inverse (effective) 
correlation length then there is a much stronger dependence of the product of Wil- 
son coefficient and matrix element of a local operator than the logarithmic one. 
Along the lines of section 3 it would therefore be important to measure the 4 quark 
correlators corresponding to the AB = 2 matrix elements of eq. (wi) in order to 
decide whether mass scales (inverse correlation lengths) occur which are dangerous 
for the heavy quark expansion. This should also be done for contributions of higher 
dimensional operators formally corresponding to higher powers in l/m&. After all, 
the inverse, effective correlation length, used in the OPE's of vacuum correlators in 
light-quark channels (set (B) in section 4), is not too small compared to m;, (10/3 
GeV <jt 4.5 GeV). 



6 Summary 

In this paper we investigated the consequences of non-perturbative coarse graining 
of operator VEV's as it was proposed in ref. ||. The focus was on light-quark 
correlators. 

After a brief review of OPE coarse graining, based on the knowledge of non- 
perturbatively calculated, gauge invariant n-point functions, we addressed the issue 
of vacuum saturation in the case of so-called disconnected diagrams. In particular, it 
was pointed out that there are ambiguities in the way approximations to the general 
prescription are implemented for coarse graining VEV's of local 4 quark operators. 
Using the simplest option for vacuum saturation, the machinery was applied to the 
light-quark correlators in the p, ax, n, and <fi channels. The spectral functions of 
the respective channels were calculated from the OPE by analytical continuation to 
time-like external momenta. Using lattice data on the gauge invariant field strength 
and scalar quark correlators, the spectra were found to be far off their experimen- 
tally measured behavior. With a 10 times smaller correlation length at dimension 



13 



6 the basic phenomenological features turned out to be contained in the spectral 
functions. However, in general the spectra are shifted to lower energies and reso- 
nance information is vague in the a\ and channels. This may be a consequence 
of excluded higher mass dimensions and the treatment of the VEV's of 4 quark 
operators. Despite the obvious shortcomings a calculation of the ratios of the first 
and zeroth moments of the spectral distributions in the a\, tt, and (f) channels gave 
quite realistic results after fixing this ratio in the p channel. Finally, we discussed 
the potential impact non-perturbative coarse graining may have on the theoretical 
determination of the inclusive, nonleptonic Ar in S s -meson decays. 
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